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Abstract 

By exact numerical and master equation approaches, we show that a central spin- 
1/2 can be configured to probe internal bath dynamics. System-bath interactions 
cause Rabi oscillations in the detector and periodic behavior of fidelity. This period 
is highly sensitive to the strength of the bath self- interactions, and can be used to 
calculate the intra-bath coupling. 
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1 Introduction 

Recent identification of potential quantum computer (QC) architectures [1] 
with long external decoherence times has shifted attention toward the negative 
effects of static internal structural flaws [2]. Isolated flawed QCs and their im- 
mediate local environments can in many cases be mapped onto self-interacting 
spin-bath models, and as a consequence spin baths are increasingly a focus of 
attention [2,3,4,5]. Studies have reported a number of interesting dynamical 
effects such as suppression of decoherence with increasing bath self-interaction 
[3,4,5,6], and large coherent shifting of subsystem dynamics [5]. 

Here we show that the coherent shift can be used to configure a single qubit 
detector to measure the strength of bath self-interactions. Structural flaws in 
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a QC arise due to variance in one-body parameters or as a consequence of 
residual two-body interactions, or even as a result of strong interactions with 
other local impurities. One body flaws can be readily detected and corrected 
in many cases. Two-body interactions are more difficult to deal with and 
their effects are difficult to predict. Knowledge of the strength of the two- 
body interactions is obviously a prerequisite for eliminating their effects. It 
has been observed that the coherent shift is sensitive to the strength of bath 
self- interaction [5]. We show that a single qubit can be configured to show 
large amplitude oscillations in fidelity when a shift is present, and that the 
period of the oscillations is highly sensitive to the strength of the bath self- 
interaction. From the period of the oscillation one can extract a quantitative 
estimate of the self-interaction strength. While we are not aware of any existing 
error correction schemes designed specifically for two-body flaws one should 
be simple to devise. Moreover, it should be possible to apply the basic ideas 
behind our detector set-up in more general contexts such as optical impurities 
in solids where knowledge of bath self- interaction could be important. 

An engaging aspect of spin-bath models is the ease with which exact results 
can be obtained for the low temperature regime relevant to quantum comput- 
ing. Unfortunately, exact simulations are necessarily limited to small baths. 
In our study this limitation is not crucial but in more general studies it might 
be very important. Consequently, we ideally would like to have a second in- 
dependent way of obtaining open system dynamics. While uncoupled spin 
baths and uncoupled oscillator baths can display similar dynamics[7], stan- 
dard boson-bath models[8,9,10,ll] are inapplicable here. This is due to the 
fact that the bath self-interactions are non-negligible in our model. We there- 
fore choose to also consider simulations obtained using a recently developed 
approximate master equation (ME) [12,13,14], which takes into account bath 
self-interactions via a mean-field approach. 

Our QC model is based on a Josephson junction (JJ) based proposal [15] which 
has an extremely long external decoherence time[15,16]. Recent experiments 
report [17,18] that the primary cause of decoherence of such JJ qubits is static 
interaction with some unidentified local two-level systems (TLSs) [19,20]. In 
one case [18] externally induced Rabi oscillations in a qubit were significantly 
distorted at times much shorter than the external decoherence time and the 
existence of the TLSs was posited as the cause. Our recent results [5] suggest 
that the distortion might be caused by a coherent shift which could arise 
just from flaws and residual interactions amoung the qubits. In any case the 
possible TLS can be viewed as distorted auxilliary qubits [19]. The computer 
and TLSs can be mapped onto a spin-bath and this is the model we address. 

We let one detector qubit interact with a statically flawed QC, and our goal 
is to show that the strength of the bath self-interaction can be measured. Our 
system is thus a single two level system and our bath consists of flawed qubits 
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and other interacting TLSs. We prepare our detector in a superposition state 
which will merely undergo phase evolution in the absence of coherent shifting. 
When shifting is present the populations and fidelity oscillate. The strength 
of the bath self-interaction is extracted from the period of these oscillations. 
We use both exact and approximate master equation simulation methods. 

While the phenomenology of decoherence and dissipation dominate the liter- 
ature of open quantum systems[3,4,5,6,7,8,9,10,ll,12,13,14] the possible ex- 
istence of coherent shifting has been acknowledged but not widely explored. 
Combinations of observables like fidelity and purity can be used to distinguish 
these coherent sorts of errors from non-unitary effects like decoherence and 
dissipation [5]. The effect of such a shift can be non-trivial and it can qualita- 
tively alter the system evolution [5]. In our detector, for example, we will see 
that pure phase evolution is changed to large amplitude Rabi oscillations. 

Distinguishing the unitary part of open system-dynamics from non-unitary 
parts is of particular importance, because unitary errors induced by the bath 
should be easily correctable. The exact Nakajima-Zwanzig ME [21] contains 
terms which result in both Markovian and non-Markovian distortions. The 
Markovian term consists of Hermitian corrections to the subsystem Hamil- 
tonian and clearly causes unitary errors. Unfortunately, this Markovian cor- 
rection cannot always be identified with the physical coherent shift. Different 
choices for the projection operator [22] lead to different Markovian correc- 
tions. It is possible in principle that a specific projection operator will have 
the property that the Markovian corrections are precisely equal to the phys- 
ical shift [12,13,14]. The identity of this projection operator is a matter of 
some interest since it is in a sense the most physical. Furthermore, the form 
of the Nakajima-Zwanzig master equation [21] which emerges from this choice 
of the projection operator would be a favorable starting point for subsequent 
approximations. 

A recently developed approximate ME [12,13,14] relies on a projection opera- 
tor for which the unitary corrections consist of a shift H$ — > Hs + BS, where 
B denotes the canonical average of the bath coupling operator B and S is the 
system coupling operator. The non-Markovian terms are all non-unitary and 
so agreement of this ME with exact calculations would strongly support this 
choice of the projection operator. We show that the predictions of the ME 
are indeed quite good and hence that this the projection operator is the most 
physical. The accuracy of the ME in this instance is of some interest since it 
also preserves positivity and has no adjustible parameters. 

Organization of this Letter is as follows. In Section 2, we explain the math- 
ematical details of our model. In Section 3 we explain our exact simulation 
method. Section 4 reviews the derivation of the exact Nakajima-Zwanzig ME 
and the approximations which lead to a solvable positivity-preserving non- 
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Markovian ME with no adjustable parameters. A numerical strategy for solv- 
ing the master equation is also reviewed. In Section 5, we present and discuss 
our exact numerical and master equation results. 



2 Self-Interacting Spin-Bath Model 



Our self-interacting spin-bath model, representing a flawed QC core consists 
of N + 1 two level systems, is given by the following total Hamiltonian 

H = H S + SB + H B (1) 



where the first term represents the system Hamiltonian 

Hs = -\b&? (2) 

the second term is a system-bath interaction operator in system S and bath 
B spaces 

SB = a^J2\^ (3) 

i=i 



and the bath Hamiltonian is given by 



1 N 

H B = --Y.(B x l ^ ) +Btaf) 
i=i 

N-l N 

+ E E (4) 

i=l j=i+l 

Note that a l a are standard Pauli spin-1/2 operators with a = {x, y, z} and the 
index i labels the qubits. 

Hamiltonian (1) is a variant of an isolated flawed QC model recently intro- 
duced to investigate the effects of one-qubit static imperfections and two-qubit 
residual interactions on the performance of QCs[2]. The actual internal dynam- 
ical effects have been explored [5] and while decoherence and dissipation are 
present, the principle errors were caused by a coherent shift. 

Here qubit will be employed as our detector while the rest represent some 
combination of flawed qubits and TLSs. Flaws will be modeled by choosing 
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the parameters of Hb randomly from some preset distribution. Note also that 
the bath qubits and TLSs interact among themselves. 

2. 1 Parameters 

The parameters used in our numerical study are based on a Josephson charge- 
qubit QC[15,16] proposal for which the experimentally accessible single qubit 
energy to perform a one-qubit rotation is Bq = le with e = 200 mK. Since 
all qubits in Hamiltonian (1) can be viewed as components of the same QC, 
the fluctuations in single qubit parameters should only differ from Bq by a 
detuning parameter, which we set to the value 8 = 0.4 e. Thus, we add noise 
to all qubits except for our detector by choosing Bf G [Bq—8/2, Bq+8/2]. We 
also make an allowance for a;- type errors by choosing Bf G [Bq—8/2, Bq + 8/2]. 
Following Ref. [2] residual two-body interactions are included by randomly 
choosing A« and uniformly from Aj G [—A, A] and G [—J x , Jx], 
respectively. While we considered a number of intra-bath coupling strengths 
J x = 0.00, 0.15, 0.50, 1.00, 2.00 in units of e, to explore the integrable to 
chaotic transition, we considered only one system-bath coupling value, which 
corresponds to the experimental value for a two-qubit rotation, A = 0.05 e. In 
all our calculations temperature is kept constant at kT = 0.25 e and time step 
of integration was At = 0.2 h/e s. 

Thus, the parameters are set so that all qubits are part of the same com- 
puter, but one qubit does not possess any one-body flaws and we use this 
as our detector. Note also that in a multi-qubit JJ QC, the existence of the 
aforementioned TLSs is expected. In this study, we will not attempt to ex- 
plicitly distinguish between qubits and TLSs because the nature of the TLSs 
is currently unknown and it is not possible to derive TLS parameters from a 
first principle. Nevertheless, the same form of pseudospin Hamiltonian used 
for the qubits is also assumed for TLSs [23] and they both have comparable 
decoherence times [23]. 

2.2 Initial Conditions 

We assume that the initial state of the composite system is of product form 
p(0)=P5(0)®p B (0) (5) 

where ps(0) = |^(0))(^(0)| is the initial pure reduced density of the system, 
and p_b(0) is the canonical bath density. We chose the system initial state as 
| -0(0)) — (|0) + |l))/\/2- In the absence of interactions with the bath this state 
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of the detector will undergo phase evolution only. But Rabi oscillations will 
appear when there are interactions with the bath. 



3 Exact Numerical Approach 

Exact numerical solutions are obtained by exploiting the low temperature 
regime for the bath degrees of freedom. Thus we approximate the initial bath 
density as 

PB(0) = T,P»\tf)(<fi\ ( 6 ) 

n=l 

where p n = exp {—E n / kT) / Y^m=i ex P {~E m / kT) are the thermal populations, 
and H B \(/>%) = E n \<j%). Note that the sum is over only the thermally popu- 
lated, lowest energy, eigenstates of the bath. Hence, n cut is a cutoff such that 
states n cut + 1 and higher are unoccupied for our given fixed temperature. 
Even though the density of bath states and thus the number of populated 
states varies slowly with J x for a fixed temperature, n cut = 20 was sufficient 
for all cases in our calculations. 

We calculate the exact reduced density ps(t) at time t via 

Tl>cut 

P s(t) = J2 P nTr B [\* n (t))(* n (t)\], (7) 

71=1 

Here the \^ n (t)) evolve according to the Schrodinger equation 

i/4|tf n (f)> (8) 



and all initial states are of the form |^ n (0)) = \tp(0))<S>\(f>n)- We used a Lanczos 
algorithm[25] for exact diagonalization of the bath Hamiltonian, for N = 10 
qubits, and an eighth order variable stepsize Runge-Kutta method [26] for the 
numerical integrations of Eq. (8). 

From the exact reduced density various measures (discussed in Section 5) can 
be calculated to determine the extent of decoherence, dissipation, and shifting. 
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4 Master Equation Approach 



Here we begin by reviewing the basic ideas behind the Nakajima-Zwanzig pro- 
jection operator technique and show how it leads to an exact but unsolveable 
master equation. Next we explain how a mean field approximation can be used 
to obtain an approximate but solveable master equation with no adjustible pa- 
rameters. 



4-1 Exact Nakajima-Zwanzig Master Equation 

Exact quantum dynamics in density operator space obeys the Liouville-von 
Neumann equation 



where the Liouville operator is L = (l/h)[Hs + SB + Hb, • }■ The Nakajima- 
Zwanzig approach begins by defining a projection operator P such that Pp(t) 
is proportional to the reduced density ps(t)- Once a P is defined, then its com- 
plement is Q = 1 — P. Operating on both sides of the Liouville-von Neumann 
equation with P and Q yields 



i^-[Pp(t)} = PLP Pp(t) + PLQ QKt) (10) 
at 

ij t [Qp{t)\=QLP Pp(t) + QLQ Qp(t). (11) 
This coupled set of equations can then be solved to find an equation for Pp(t): 



i^-P P (t)=PLP Pp(t) + PLQ e - m & Qp(0) 
at 

t 

-1 J dt'PLQ e «3i0(f-t) QLP Pp{t') (12) 


Tracing over the bath degrees of freedom at this stage would lead to a general 
master equation for the reduced density. 

The properties of the resulting master equation will in part be determined by 
the choice of P. In particular, it may be possible to choose a P such that the 
non-Markovian term generates non-unitary evolutions only. Assuming that the 
system and thermal bath are initially uncorrected (i.e. p(0) = ps(0) ® Ab(0)) 
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and that the number of bath degrees of freedom is large enough so that the 
system does not significantly alter the equilibrium state of the bath, then the 
total density pit) at time t will be approximately given by ps(t) ®Ab(0). Thus, 
it seems to make sense to define the projection operator P via 

Pp(t) = p s (t) ® p B (0). (13) 



It should be noted that this was not the actual projection operator employed 
by Zwanzig[21]. This particular definition has however been widely employed 
in the literature, and it is well known to be of non-Hermitian type[12, 13, 14,22]. 
With this P the master equation reduces to the form 

±p s (t) = -t[H s + BS,p s (t)} (14) 
t 

-J dt' Tr B {p B (0) LQe 1 ®^^' _€) QL}p s (t 1 ) 
o 

in which we see the emergence of the coherent shift Hs — > Hs + BS. The 
problem with this master equation is that evaluation of the memory kernel is 
extremely difficult. 



4-2 Non-Markovian Mean Field Master Equation 



More useful, although approximate, equations are obtained by making approx- 
imations for the memory kernel. Here we explore a mean field approach. 

Since P and Q are not Hermitian the operator QLQ, mediating correlations 
between system and bath, is not Hermitian either [13,14]. Assuming that QLQ 
has a complete complex spectral decomposition, then for evolution forward in 
time t > we write 

qz4 = £K-^)W)WI ( 15 ) 



where |0+) and (<&t\ are the right and left eigenvectors, respectively. Substi- 
tuting (15) back into (12) and tracing over the bath gives 



d . , 



h 



H s + BS,p s (t)] - dt'JC(t-t') Ps (t') 



(16) 



where 
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t{t-t')= ^ e -»(^-%)(w') 

xTr B {p B (O)LQ|0+)($+|QL}. (17) 

Further progress now requires some form of approximation for above sum 
since the eigenspectrum cannot be obtained in practice. A mean field approach 
seems logical if |0+) and (<&+| are random functions as they would be for a 
chaotic bath. Further assuming that the eigenvalues — 27^ are statistically 
independent of the eigenvectors then suggests that the sum should be replaced 
by its average 

(e-^-^)*> ]T |0+)($+| = (cosM)e^)i 5 ® i B . 
Defining W^(t) = (cos(cjt)e~ 7 *) then we obtain 



j t Ps(t) = ~[H S + 5S,p 5 (i)] - J dt'W(t - t')t DPs {t'). (18) 



t 



Here Cd = (C/h 2 ){[psS, S] + [S, Sps)]} is a dissipative Lindblad-Kossakowski 
superoperator [24], and and C = B 2 — B 2 is the variance of the bath coupling 
operator. Here the non-Markovian term is purely dissipative and the only 
unitary effects arise from the shift BS. 

In the Markovian limit (i.e. when memory effects are negligible) 



l 00 
J dt'Wit - t')C D p s (t') ~ J dt'W(t')C D p s (t) 



~T B £ DPs {t) (19) 

where tb is the (short) memory relaxation time. In this limit Eq. (18) reduces 
to the standard completely-positive-dynamical-semigroup (CPDS) form[24] 

j t Ps{t) = - l -[H s + SB,p s (t)\ - r B C D p s (t). (20) 

We however will employ the non-Markovian form since memory effects are 
important in our system. Analysis of the spectral properties of QLQ suggests 
that an appropriate memory function for chaotic baths is of the form[14] 



W{t) = [1 - A(pt) 1 + l -{ptf - J-(pt)3] e -^) 2 /8 (21) 



9 



where expressions for the quantities p and q are presented elsewhere [27]. Here 
W(t) is positive and real. 

It is important to note that the master equation (18) preserves positivity 
[12,13,14] for our positive memory function, and under certain conditions out- 
lined in [28] it also preserves complete positivity. The master equation is free 
of adjustable parameters and thus it can be solved numerically without any 
further approximation. Hence the master equation (18) can be viewed as a 
non-Markovian generalization of CPDS theory [24], which draws a sharp dis- 
tinction between unitary and non-unitary contributions to the system dy- 
namics. New types of projection operators, leading to other non-Markovian 
generalizations of CPDS theory [24] are discussed in [29]. Agreement of our 
master equation with exact calculations would then suggest that Eq. (13) is 
indeed the physically meaningful choice for P. 



4-3 Numerical Solution of Master Equation 



We use a recently developed numerical method for solution of non-Markovian 
equations [30]. The method is implemented by converting integro-differential 
equations to ordinary differential equations. First, define a space-like time 
variable u and a smoothed density function 

t 

X(t, u) = f(u) J dt' W{t -t' + u)p(t'), (22) 
o 



where f(u) is a damping function with /(0) = 1. Then, substituting Eq. (22) 
into Eq. (18) we obtain two coupled ordinary differential equations for p(t) 
and x(t,u) 



f t Ps(t) = -^H s + SB,p s (t)] 

-^{[x(t,0)S,S] + [S,Sx(t,0)]}, (23) 

±x(t,u)=f(u)W(u)p s (t) 

dx(t,u) f'(u) . /csA , 
du f(u) 

which can be solved using standard methods. 

For the spin-spin-bath model of our study these equations have the form 
(h=l) 
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-2C{ X (t,0)-a^ X (t,0)a x ^} (25) 
j t x(t,u)=e-° u2 W(u)p(t) 

+ ^^ + 2gux(t,u), (26) 

where B = E x and C x = Y? x — Y? x . Here E x = X)fcLi and the overbar 

denotes a canonical average over the initial bath state. The parameters of 
the memory function, C x , and E x are calculated using exact energies and 
eigenvectors of the bath Hamiltonian. 

We found that a damping function of the form f(u) = e" 9 " 2 with g = 
9.9/[(n — I) At} 2 , as suggested in Ref. [30], was accurate for our problem. Using 
W(u) = W(\u\) for negative values of u the auxiliary differential equation 
(26) is solved on a grid of points uj = (n + / — j) At with j = 1, . . . , n and 
I = mt(.338n). Here At = 0.2 h/e is the time-step employed in the dynam- 
ics. Partial derivatives with respect to u are calculated by using a discrete- 
variable[31] matrix representation. Converged results were obtained for n = 40 
grid points. The ordinary differential equations (25) and (26) were integrated 
by an eighth order variable stepsize Runge-Kutta routine[26]. 



5 Results 



We use purity and fidelity to detect and quantify deviations from pure phase 
evolution in our detector. Purity, defined by V(t) = Tr s[ps{t)] 2 , is a good 
measure of decoherence and dissipation since it is insensitive to the coherent 
shift. For pure initial states, such as we will employ, the dynamics of the system 
will ideally have a purity equal to 1. This would be the case if interaction with 
the bath causes only coherent shifting but no non-unitary effects. 

Fidelity, defined by F(t) = Tr s [p s (t)p l g eal (t)}, is sensitive to both unitary 
and non- unitary effects. Here p l g eal (t) is the system density at time t in the 
absence of system-bath interactions. Fidelity measures the proximity of an 
open system's evolution to its free evolution. Hence, for pure initial states the 
ideal value of fidelity is also 1. Comparison of the magnitudes of purity and 
fidelity gives an indictation of the presence of coherent shifts. 

Recall now that we have chosen a initial state (|0) + 1 1))/ \[2 for the system. In 
the absence of coherent shifting this state will undergo only phase evolution 
without significant population transfer. In the presence of coherent shifting 
we expect that the system will undergo large amplitude Rabi oscillations. We 
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Fig. 1. (Colour online) Exact numerical results are plotted in (a) for V(t) and in 
(b) for F{t). Master equation results are plotted in (c) for V{t) and in (d) for ^{t). 
Five different intra-bath couplings J x = 0.00, 0.15, 0.50, 1.00, 2.00 are shown in each 
subplot. 

thus also look at the populations of the states to additionally monitor the 
emergence of the coherent shift process. 



5.1 Short Time Dynamics 

The short time behaviors of V(t) and J-"{t) are seen in Figure 1 for both exact 
and master equation results for five different intra-bath couplings, J x . The 
systematic improvement in V(t) with increasing J x is observed for both ex- 
act and master equation calculations. For the integrable bath i.e. J x = 0.00 
the decoherence is at its maximum. Above J x = 0.15 chaos sets in and the 
decoherence is systematically reduced. For strong chaos with J x = 2.00 the 
decoherence is almost totally suppressed. This result is in agreement with ear- 
lier studies[3,4,5,6] in which the bath chaos is predicted to reduce decoherence 
and dissipation. 

Partial recurrences of the purity are observed in the exact calculation in the 
integrable regime. This is obviously a consequence of the small bath size and 
the effect is not observed in the master equation results. Agreement between 
the master equation and exact calculation clearly improves with increasing 
J x , and this is a consequence of the fact that the dynamics becomes more 
Markovian with greater bath chaos. Markovian character requires a separation 
of relaxation times for the system and bath degrees of freedom [32] . A chaotic 
bath can relax internally, but the only available relaxation mechanism for an 
integrable bath is through the system. Thus, chaotic baths tend to be more 
Markovian. Since we assumed the bath was chaotic in its derivation we should 
not be surprised that the ME is more accurate for larger J x . 
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(a) 




Fig. 2. Canonical average (a) and variance (b) of bath coupling operator plotted vs. 

Jx- 

The dissipation term in Eq. (18) is preceded by a factor which is equal to the 
canonical variance in the bath coupling operator. This variance is a measure of 
the number of bath states which can contribute to decoherence and dissipation 
of the system. Fig. 2(b) shows that this variance declines with J x and so this 
in part explains why Fig. 1 shows an improvement in purity with increasing 
J x . It has been argued elsewhere [3] that the off-diagonal matrix elements of 
the bath coupling operator tend to vanish as the chaotic regime is approached. 
This lack of selection rules is largely responsible for the decline in the variance. 

The fidelity plots of Fig. 1(b) show the presence of much larger deviations from 
unity than that seen in the purity plots. This supports the conclusion that 
large coherent shifts are again occuring as they did in our previous study [5]. 
This shift is clearly much more harmful than decoherence or dissipation as 
an potential error source for quantum computation [5]. We will see that the 
coherent shifting can be put to a good purpose, however. It is also worth noting 
that while the coherent shift vanishes in some standard models of decoherence, 
it is nonzero in general for self-interacting baths. Thus, we expect that shifts 
of this type and magnitude would also arise for solid state baths. 

The fidelity plots show a quite different trend with increasing J x . After a short 
time F{t) begins to decay more rapidly for larger values of J x . This effect is 
well captured by the master equation except, again, for the smallest J x values. 
In Fig. 2 (a) we plot the canonical average of the bath coupling operator. It 
clearly grows in magnitude with J x and so the shift SB, predicted by Eq. (18), 
should also grow. This is clearly in agreement with the trend in Fig. 1(b). This 
sensitivity of fidelity to J x also suggests that the shift is sensitive to the bath 
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Fig. 3. (Colour online) Matrix elements of system density operator ps(t) for 
J x = 1.00: Exact numerical results are given by solid lines and master equation 
results by dots, (a) Diagonals of ps(t), (0\ps(t)\0) (blue) and (l\ps(t)\l) (black), 
(b) Real part of off-diagonal element of ps(t), Re{ (Ol^ (t) 1 1)} (green) and imagi- 
nary part of off-diagonal of ps(t), lm{(0\ps(t)\l}} (red). 

self-interaction strength. 

In Figure 3 we plot the matrix elements of the detector qubit density ps(t) 
for J x = 1.00. Recall that only phase evolution is expected in the absence 
of coherent shifting. Only very small population transfers would be expected 
as a result of dissipation. In Fig. 2(a) we see large Rabi oscillations in the 
populations which are thus a direct consequence of coherent shifting. In ad- 
dition, complete agreement between the master equation and exact numerical 
results is evident for populations (see 3(a)) and phases (see 3(b)). This evi- 
dence strongly suggests that our choice (13) for the projection operator P is 
the physically correct one. 



5.2 Very Long Time Dynamics 

In Figures 4(a) and 4(b) we plot exact numerical results for purity and fidelity, 
respectively, for five different values of J x . The master equation results for the 
same quantities are plotted in Figures 5(a) and 5(b). 

Figures 4(a) and 5(a) show that suppression of decoherence with increasing J x 
also occurs in the long time dynamics of the purity. In addition, the purity plots 
for the exact calculations reflect the presence of long time partial recurrences 
in the integrable bath regime. These recurrences are almost certainly caused 
by memory effects induced because of the small bath size and they do not 
appear in the master equation results. The magnitudes of the exact and ME 
predictions for the purity are however in good agreement. With increasing 
bath chaos the recurrences are considerably reduced and agreement between 
the exact and ME results is much better. 
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Fig. 5. (Colour online) Master equation results for long time dynamics: (a) Purity 
V(t) vs. time, (b) Fidelity F(t) vs. time for J x = 0.00,0.15,0.50, 1.00,2.00. 

In Figures 4(b) and 5(b) we plot the exact numerical results and master equa- 
tion results, respectively, for the long time dynamics of the fidelity for five 
different values of J x . In the long time limit, the unitary effects of system- 
environment interactions overwhelm the non-unitary ones, and the contribu- 
tions of decoherence and dissipation to the dynamics are hardly noticeable 
anymore. This is true even in the integrable regime. Small magnitude, high 
frequency, oscillations are still noticeable. However, the fidelity plot displays 
an additional long time, large amplitude, periodicity. Moreover, the period of 
the fidelity is strongly dependent on the magnitude of J x . There is again a 
discrepancy between the ME and exact results in the integrable regime, but 
the agreement is excellent in the chaotic cases. 

Explaining the sensitivity of the fidelity period to J x is quite straightfor- 
ward. Neglecting the effects of decoherence and dissipation it can readily be 
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argued that the shifted system dynamics should beat with a frequency of 
f2 = \J B 2 + Bq 2 /A/h while the unperturbed system phase evolves with fre- 
quency iv = Bq/2H. The period of Rabi oscillations in Fig. (3) is ir/fl. It then 
follows that the fidelity takes the form 



The second term is very small since B 2 is very much smaller than f2 2 . The 
third term is small since Q — B^/2 < Q, and the fourth term is of order 1. 
Hence, the small magnitude oscillations in fidelity have frequency 2(Q + a;). 
The large amplitude oscillations are caused by the fourth term, and they have 
period 7r/ (Q — uo) ~ }iBq/2B 2 . Generally, B is much smaller than Bq but B 
increases with J x resulting in a shorter period. This gives rise to the changes 
in the period of the fidelity and fully explains the behavior observed in Figs. 
4(b) and 5(b). Since 2(Q + u) varies more slowly with B the Rabi oscillations 
are not very sensitive to J x . 

Thus, by measuring the period of the fidelity oscillation one can obtain an 
estimate of B, and from Fig. 2(a) we can then obtain the magnitude of J x . 
Hence, we have a detector of the strength of bath self-interaction for this JJ 
QC model. The same basic setup should also carry over to the case of oscillator 
baths, and the technique could potentially be used to measure the strength of 
anharmonic interactions in the solid state. 



6 Summary 

We studied the decoherence dynamics of a central spin interacting with a self- 
interacting spin bath by exact numerical calculations and a master equation 
approach. We obtained good agreement between the two methods. The model 
represents a qubit detector interacting with an isolated flawed QC. In the ab- 
sence of interactions with the QC the detector undergoes phase evolution only. 
When the detector experiences a coherent shift as a result of interaction with 
the QC it begins Rabi oscillations. The fidelity also exhibits a periodicity on a 
much longer timescale with a period which is sensitive to the strength of bath 
self-interaction. Measuring the period allows one to find the intra-bath cou- 
pling strength. The same approach could be used for measuring anharmonic 
interactions in more general contexts. 




cos2(f2 + iv)t 



(27) 



16 



Acknowledgements 



The authors acknowledge the support of the Natural Sciences and Engineering 
Research Council of Canada and the computer resources provided by West- 
Grid. 



References 



[1] M.A. Nielsen, I.L. Chuang, Quantum Computation and Quantum Information, 
Cambridge University Press, Cambridge, 2000. 

[2] B. Georgeot, D.L. Shepelyansky, Phys. Rev. E 62 (2000) 6366; B. Georgeot, D.L. 
Shepelyansky, ibid. 62 (2000) 3504; G. Benenti, G. Casati, D.L. Shepelyansky, 
Eur. Phys. J. D 17 (2001) 265. 

[3] L. Tessieri and J. Wilkie, J. Phys. A 36 (2003) 12305. 

[4] CM. Dawson, A.P. Hines, R.H. McKenzie, G.J. Milburn, Phys. Rev. A 71 
(2005) 052321; S. Paganelli, F. de Pasquale, S.M. Giampaolo, Phys. Rev. A 
66 (2002) 052317; M. Lucamarini, S. Paganelli, S. Mancini, Phys. Rev. A 69 
(2004) 062308; R. Alicki, Open Sys. & Inf. Dyn. 11 (2004) 53; X.S. Ma, A.M. 
Wang, X.D. Yang, H. You, J. Phys. A 38 (2005) 2761; X.-Z. Yuan, K.-D. Zhu, 
Europhys. Lett. 69 (2005) 868; X.-Z. Yuan, K.-D. Zhu, Z.-J. Wu, Eur. Phys. 
J. D 33 (2005) 129; T. Prosen, T.H. Seligman, J. Phys. A. 35 (2002) 4707; M. 
Znidaric, T. Prosen, J.Phys. A 36 (2003) 2463. 

[5] M. Cetinba§ and J. Wilkie, submitted for publication to Phys. Lett. A. 

[6] A.S. Sanz Ortiz, Y. Elran, P. Brumer, J. Wilkie, and M. Getinba§, manuscript 
in preparation. 

[7] A. A. Golosov, S.I. Tsonchev, P. Pechukas and R.A. Friesner, J. Chem. Phys. 
Ill, 9918 (1999). 

[8] CW. Gardiner and P. Zoller, Quantum Noise, 3rd Ed., Springer, Berlin, 2004. 

[9] C. Anastopoulos, B.L. Hu Phys.Rev. A 62 (2000) 033821; A.O. Caldeira, A.J. 
Leggett, Ann. Phys. 149 (1983) 374; R.P. Feynman, F.L. Vernon, Jr., Ann. 
Phys. 24 (1963) 118; F. Haake, R. Reibold, Phys. Rev. A 32 (1985) 2462; A.O. 
Caldeira, A.J. Leggett, Physica A 121 (1983) 587. 

[10] L. Diosi, Quantum Semiclass. Opt. 8, 309 (1996); W.T. Strunz, Phys. Lett. A 
224, 25 (1996); L. Diosi and W.T. Strunz, Phys. Lett. A 235, 569 (1997); L. 
Diosi, N. Gisin and W.T. Strunz, Phys. Rev. A 58, 1699 (1998); W.T. Strunz, 
L. Diosi and N. Gisin, Phys. Rev. Lett. 82, 1801 (1999); J. Gambetta and HM. 
Wiseman, Phys. Rev. A 68, 062104 (2003). 



17 



[11] H.-P. Breuer, F. Petruccione, The Theory of Open Quantum Systems, Oxford 
University Press, New York, 2003. 

[12] J. Wilkie, Phys. Rev. E 62 (2000) 8808. 

[13] J. Wilkie, J. Chem. Phys. 114 (2001) 7736. 

[14] J. Wilkie, J. Chem. Phys. 115 (2001) 10335. 

[15] J. Q. You, J.S. Tsai. and F. Nori , Phys. Rev. Lett. 89 (2002) 197902. 

[16] Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 73 (2001) 357 and 
references therein. 

[17] Y. Nakamura at al Nature (London) 398 (1999) 786; J.M. Martinis at al, Phys. 
Rev. Lett. 89 (2002) 117901; D. Vion at al, Science 296 (2002) 886; Y. Yu at 
al, Science 296 (889) 2002; I. Chiorescu at al, Science 299 (2003) 1869. 

[18] R.W. Simmonds, K.M. Lang D.A. Hite, S. Nam, D.P. Pappas, and J.M. Martinis 
Phys. Rev. Lett. 93 (2004) 077003. 

[19] S. Ashhab, J.R. Johansson, and F. Nori, Physica C 444 (2006) 45; A. Shnirman, 
G. Sch6n, I. Martin, and Y. Makhlin, Phys. Rev. Lett. 94 (2005)127002; L.-C. 
Ku and C.C. Yu, Phys. Rev. B 72 (2005) 024526. 

[20] S. Ashhab, J.R. Johansson, and F. Nori, New J. Phys. 8 103 (2006); Y.M. 
Galperin, D.V. Shantsev, J. Bergli and B.L. Altshuler, Europhys. Lett. 71 (2005) 
21; F. Meier and D. Loss, Phys. Rev. B 71 (2005) 094519. 

[21] S. Nakajima, Prog. Theor. Phys. 20 (1958) 948; R. Zwanzig, J. Chem. Phys. 33 
(1960) 1338; R. Zwanzig, in Lectures in Theoretical Physics, Vol. 3, Interscience, 
New York, 1961. 

[22] A. Suarez, R. Silbey, I. Oppenheim, J. Chem. Phys. Ill (1992) 5101; C. Meier, 

D. J. Tannor, J. Chem. Phys. Ill (1999) 3365; V. Romero-Rochin, A. Orsky, I. 
Oppenheim, Physica A 156 (1989) 244; W. Peier, Physica (Utrecht) 57 (1972) 
565. 

[23] A.M. Zagoskin, S. Ashhab, J.R. Johansson, and F. Nori, Phys. Rev. Lett. 97 
(2006) 077001. 

[24] G. Lindblad, Commun. Math. Phys. 48 (1976) 119; V. Gorini, A. Kossakowski, 

E. C.G. Sudarshan, J. Math. Phys. 17 (1976) 821; R. Alicki, K. Lendi, Quantum 
Dynamical Semigroups and Applications, Springer, Berlin, 1987. 

[25] R.B. Lehoucq, D.C. Sorensen, C. Yang, ARPACK Users' Guide: Solution of 
Large-Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods, 
SIAM, Philadelphia, 1998. 

[26] E. Hairer, S.P. Norsett, G. Wanner, Solving Ordinary Differential Equations 
I. Nonstiff Problems, 2nd Ed., Springer Series in Computational Mathematics, 
Vol. 8, Springer- Verlag, Berlin 1993. 

[27] M. Cetinba§ and J. Wilkie, manuscripts in preparation. 



18 



[28] A. A. Budini, Phys. Rev. A 69 (2004) 042107. 

[29] H.-P. Breuer, Phys. Rev. A 75 (2007) 022103; H.-P. Breuer, J. Gemmer, and 
M. Michel, Phys. Rev. E 73 (2006) 016139. 

[30] J. Wilkie, Phys. Rev E 68 (2003) 027701. 

[31] D.T. Colbert, W.H. Miller, J. Chem. Phys. 96 (1992) 1982. 

[32] Y.C. Cheng and R.J. Silbey, J. Phys. Chem. B 109 (2005) 21399. 



19 



